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Abstract: We examine methods that have been proposed for determining the he- 
licity structure of decays of new resonances to third generation quarks and/or lep- 
tons. We present analytical and semi-analytical predictions and assess the appli- 
cability of the relevant variables in realistic reconstruction scenarios using Monte 
Carlo-generated events, including the effects of QCD radiation and multiple parton 
interactions, combinatoric ambiguities and fast detector simulation. 
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1. Introduction 

The third generation is thought to be intimately connected to the mechanism respon- 
sible for electroweak symmetry breaking due to the comparatively large top, tau and 
bottom Yukawa couplings. New, third generation resonances are further favoured in 
comparison to those of the first and second generation by flavour constraints, coming 
from experiments investigating flavour- changing processes [1]. 
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Collider searches for new resonances decaying to third generation quarks or lep- 
tons (or both in the case of leptoquarks) are particularly challenging. If the new 
resonance decays leptonically to either a r or u T , the final state necessarily pos- 
sesses missing energy, which prevents us from directly reconstructing the events. If 
it decays to b quarks, the signal has to compete with the high QCD cross sections, 
which can imitate the topology. In the case that the new resonance decays to top 
quarks, these may subsequently decay to several jets (for the hadronic decay), or to 
a jet plus lepton and missing energy (for the semi-leptonic decay). Both top quark 
decay modes require further reconstruction, where combinatorial ambiguities, QCD 
backgrounds and the unknown neutrino momenta can be problematic. 

Despite the aforementioned difficulties, third generation resonances provide us 
with an interesting opportunity. The decays of the resulting third generation fermions 
(apart from u T ) may allow us to determine the chirality structure of their couplings 
to the new resonance. To fully determine the underlying theory and reconstruct the 
Lagrangian terms, determination of the coupling structure, as well as the spin of the 
new resonance, is necessary. The determination of the helicity of top quarks and r 
leptons is made possible by the fact that the angular distributions of their daughter 
particles are highly correlated to the helicity of the parents. Determination of the 
helicity of top quarks has been investigated theoretically in detail [2-4]. Several 
variables have been proposed and QCD corrections to these have been calculated. 
Similar variables have also been proposed for r leptons [5-7]. The b quark hadronizes 
before decay and produces B mesons which have relatively long lifetimes and thus 
produce displaced vertices. This allows for tagging those jets that originate from b 
quarks but washes away the effect of its helicity from the angular distributions of the 
associated daughter particles. Thus, the helicity couplings of resonances containing 
a b quark has to be inferred by first determining the spin of the parent and sister 
particles. 

In this paper we investigate the feasibility of helicity measurement at the Large 
Hadron Collider (LHC) for top quarks and r leptons resulting from the decays of new 
third generation resonances. We first review the variables that have been previously 
defined for top and r decays and reproduce the relevant distributions, comparing 
them to results from a general-purpose Monte Carlo event generator, HERWIG++ [8] 
(section In the case of hadronic top quark decays, the variables have thus far 
mostly been considered in the highly-boosted case. We relax this approximation and 
attempt to determine their usefulness in more realistic reconstruction situations. For 



simplicity, in sections |3J] and [O] respectively, we examine the applicability of these 
methods to a model containing a new heavy vector boson, a Z 1 ', which possesses 
decays to the third generation: either to a top and a light quark (specifically tu, 



1 ln appendix |A] we also present a set of angular variables that complement the energy fraction 
variables given in the literature thus far. 
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iu) 2 or a pair of r leptons. We consider reconstruction of these topologies using 
mass-shell constraints which lead to polynomial equations. In the latter case we 
also employ the r lepton decay vertex information. Subsequently, in section [373|, we 
consider the more challenging case of a third-generation leptoquark model, focusing 
on pair-production of these followed by decays to a top quark and a r lepton. The 
reconstruction technique we employ is related to those presented in Ref. [13]. Finally, 
we present our conclusions in section 0L 



2. Variable definitions 

In the decays of top quarks and r leptons, the distribution of the cosine of the angle 
between the spin axis of the parent particle k and the daughter particle i, cos(9j, is 
given by: 

P(cos0<) = ^(1 + P kKi cos 6i) , (2.1) 

where Pk = ±1 corresponds to the spin-up or spin-down states of the parent respec- 
tively and Ki is known as the 'resolving power' of i. The decay product i can also be 
taken to be a jet, i.e. containing more than one particle. The resolving power has 
been calculated in higher orders and for various decay products [14-17]. 

2.1 Daughter-to-parent energy ratios, x Pi i 

If tops and taus are produced as a result of the decay of a heavy resonance, their 
own decay products may be collimated due to their large velocities. In this case, it 
is useful to define the energy fraction x^i = £k,i/£k between the parent particle k 
(= t, t) and the one of its daughter particles, i, rather than the angle between them. 
These are formed in the laboratory frame for both the tau leptons and top quarks 
as: 

•^T,jet £')ct/ &t i 

^top,b = £b/£top , (2.2) 

respectively. For brevity, we will write x top and x T to denote the preceding variables. 

Analytic predictions for these variables can be derived using Eq. ( |2.1|) , by trans- 
forming from the angular variable cos 9i (defined in the centre-of-mass frame) to the 
energy ratio (defined in the lab frame). In the highly-boosted cases, where the boost 
factor, (3k = \p k \/Ek, is taken to be unity, their distributions are shown in Fig. [I] for 
the r lepton decay mode r — > ixv T (left) and hadronic top (right). The distribution 
of x top red has a cut-off at the maximum value of x t0 p,max = 1 ~ m vt-V m to P ~ 0-79 
and that of x T at a minimum value of x Tjm ; n = m^/m 2 ., as a result of the kinematic 

2 Note that such models have been proposed as possible explanations to the Tevatron tt asym- 
metry [9-12]. 
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restrictions imposed by the mass of the W and the pion respectively. In the case of 
the r, we show only the r — > tiv t mode, for which m n ~ 0.14 GeV, which results in 

a Small ^7-jmiri' 

The exact analytic forms of the distributions are given by [4]: 3 

1_ diV = 1 m 2 T ( _ 1 P T (m 2 T + m 2 ct ) + 1 2P T m 2 T \ 
N dx T /3 T m 2 - m 2 et \ /3 T m 2 T - m 2 ct (5 T m 2 - m 2 ct T J 
1 dN 1 m 2 ( 1 1 2m 2 \ 

TV dxtop Pt mf op - m z w \ (3 t p t m^ op - mfo J 

where Pi = ±1 (for % = r,t) represent right or left helicities of the r or top and: 

Kb ~ m 2 op + 2m 2 w - °- 4 ' (2 ' 5) 

is [4] the resolving power of the 6-quark at leading order, 4 m top , rnw and m T are the 
top quark, W boson and r lepton masses respectively. The approximation — > 1 is 
almost always good for r leptons, and with the current beyond-the-standard model 
(BSM) third generation limits rising as the LHC experiments produce more exclusion 
regions (see, e.g. [19,20]), it should be even more applicable for most BSM models 
that include heavy particles that decay to tops. Note that in the opposite limit, 
of (3 t — > 0, the above distributions tend to delta-functions, S(xk — Xk,o), since the 
lab frame and parent top (or tau) decay frames become identical. The values of 
Xkfl, determined by the masses of the decay products in the two-body decay, are 
x Tj0 — (m 2 . — ml)/(2m T ) and x t o P ,o — ( m to P — m w) I \2m tov ) for the tau and top 
distributions respectively, neglecting jet and bottom quark masses. 

It is obvious from Eqs. ( [2.3| ) and ( |2.4| ) that the effect of j3 t < 1 is to alter the 
predictions for the left- and right-handed distributions. The result of this effect is 
shown in Fig. |2], where we show the x top right-handed helicity distributions for fixed 
Pt = 1-0, 0.9, 0.8, 0.5 and 0.1. It is obvious that for /3 = 0.9 and = 0.8, the variation 
from the (3 t = 1 line is small, whereas for j3 t — 0.5 and f3 t = 0.1, the slope and position 
of the lines is dramatically altered. 5 This is due to the fact that the slope of Eq. (|2.1| ) 
is inversely proportional to fi 2 . One may then safely consider the region (3 t > 0.9 as 
the 'highly-boosted' region in the case of top quark decay. This would correspond to 
top quarks with energies E t > 2.5m top . In realistic scenarios in hadron colliders, (3 t 
is not fixed, and one would need to integrate over the allowed values to obtain the 
full result. 



3 Note that in [4], Eq. ( 2^4 ) is missing a factor of ^ — in the second term. The corresponding 
distributions, however, appear to have been constructed with the correct formulae. 

4 The authors of Ref. [18] present a study of polarisation observables at next-to-leading order 
with parton showers in H~t and W~t production. They find that in those cases these observables 
are robust against higher-order corrections. 

5 Note that there is also a lower limit on allowed values of Xtop due to kinematics. The limits for 
arbitrary top quark boosts /3 are given by (1 - (3 t ) < m x top < (1 + Pt)- 

'"top " l W 
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Figure 1: The predictions for the energy fractions x T = £- jet /£ T for the decay r — > ttu 
(left) and xt op = £b/£top (right) in the highly-boosted cases (fi T = 1 and j3 t = 1). 



20 



5 



1 ' 1 1 ' 1 1 1 1 1 1 ' 1 1 1 ' 1 1 

\ 

\ 

— 1 tR, fi 


_ 


tR, /S 


= 0.9 - 


\ ~ ■ ~ tR, 


= 0.8 - 


I tR, j8 


= 0.5 - 


' tR.fi 

I 


= 0.1 - 











0.0 0.2 0.4 0.6 0.8 1.0 



Figure 2: The predictions for the energy fraction x top = £b/£top for the cases fit = 
1.0,0.9,0.8,0.5 and 0.1. 



The r lepton can decay via multiple channels and each of these channels con- 
tributes to the total distribution. We do not attempt here to reproduce its analytic 
form, as this would require calculating the distributions of Eq. (|2.3| ) corresponding to 
each decay mode and integrating over the distribution of the mass of the r jet, m- ]et) 
with certain weight obtained by the matrix element squared. Instead we present 
results constructed from the Monte Carlo-simulated decays of the r lepton using the 
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Figure 3: The Monte Carlo predictions for the energy fractions x T = Ei e t/S T for all the 
decay modes when the r is highly boosted (f3 T = 1). 



HERWIG++ event generator. We show the resulting distributions for the left- and 
right-handed highly-boosted taus in Fig. ||. The contributions of the multiple decay 
modes clearly modify the behaviour compared to that of the single decay mode which 
appears in Fig. |]. 

2.2 Semi-leptonic top variable, u 

In the case of semi-leptonic top decays, t — > biug, where a fe-jet is tagged and an 
electron or muon is identified, we can calculate the fraction of the visible energy 
carried away by the lepton, u: 

u = tA^t (2.6) 

The resulting distributions of the variable u for highly-boosted — > 1) [4] left- and 
right-handed tops are shown in Fig. |j. Highly-boosted Monte Carlo-generated curves 
are also shown for comparison, with and without final state radiation (FSR). 6 The 
kink at u = m^/m^ op ~ 0.215 is due to the fact that there exists a minimum possible 
value of the lepton energy in the top rest frame, given by E ijrain = m^/(2m top ), which 
arises when the lepton is anti-aligned with the top boost direction. The maximum 
value of the energy is -E^^ax = "^top/2 and arises when the lepton is aligned with the 
boost direction. This is clarified in Fig. j|, where the schematic diagram demonstrates 
the decay of a top in its rest frame. 



6 Note that the curves of Fig. 2 in Ref. [4] seem to fit the FSR on case rather than the one 
obtained by integrating the differential width directly, which should fit the FSR off case. 
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Figure 4: The distribution 1/rdr/du of the fraction of visible lab frame energy carried 
by the lepton in a highly-boosted semi-leptonic top (i.e. (3 t = I), u = (££ + £{,) is shown. 
The blue curve and red curves represent left- and right-handed top quarks respectively. 
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Figure 5: A schematic diagram of the decay of the top in its rest frame into a W and a 
6-quark, with subsequent decay of the W to a lepton and neutrino. The two configurations 
shown correspond to the minimum and maximum energy configurations of the lepton, 
corresponding to E^^ min = m^/(2m t0 p) and Eg >max = m top /2 respectively. The minimum 
energy of the lepton causes the kink in the u variable distribution. 

It is important to note that the variable u has the advantage that there is no need 
to explicitly reconstruct the top quarks in order to form it, even in the case of /3 t ^ 1. 
As a result, it is expected to be less sensitive to the reconstruction systematics that 
may enter other energy fraction variables. 



3. Applications 



We examine a model of a heavy vector boson (Z') decaying to tops or taus and 
a scenario of pair-production of third generation scalar leptoquark states. In what 
follows we will neglect neutrino masses, and use the following values for the top quark, 
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W boson and r lepton masses respectively: m top = 174.2 GeV, mw = 80.403 GeV 
and m T = 1.777 GeV. 



3.1 Flavour-changing Z' 

We first examine the application of the variables x top and u that have been defined 
in the previous sections, on a model of a Z' boson that possesses flavour- changing 
quark couplings, described by the Lagrangian density: 

C = gf, R Z^u RllM t R + g% L Z'»u Ll ,t L 

+ g^ R Z'^u Rl ,u R + g u £ L Z'»u Llil u L + h.c. , (3.1) 

where g 1 ^ and g l ^ R are the flavour-changing left- and right-handed parameters that 
we will be varying. 

It is evident that the Z' bosons of this model would be produced at a hadron 
collider via light uu initial states. They will decay to both light uu and tu (and ut). 
After discovery of such a state, for example in dijet resonance searches, determining 
the helicity structure of the uuZ vertex will be challenging, if not impossible. Hence 
one would concentrate on determining the helicity structure of the vertex which 
involves the top quark. 

3.1.1 Parton-level results 

We consider a Z' described by the particular model given by the Lagrangian of 
Eq. ( ft.ip , of mass 1.5 TeV, choosing either a purely left-handed third generation 
coupling: g^i R = 0, g 1 ^ = 1 or a purely right-handed one: g l % R = 1, g t ^ l , L = 0, 
keeping in both cases g^ R = g^f L = l. 7 We show Monte Carlo results of the 
distributions of the variables x top and u, for a 14 TeV LHC, in the dashed curves in 
Fig. H These results were obtained using parton-level events, ignoring initial- and 
final-state radiation, hadronization, and applying no rapidity or momentum cuts on 
the particles. The plots in Fig. |6| also contain semi-analytic predictions that take 
into account the finite boost of the top quark in the lab frame (solid curves). The 
semi-analytic x top distribution was produced by assuming that the cosine of the 
angle of the emitted 5-quark in the top quark rest frame, cos 6^, was distributed 
according to Eq. ( |2.1| ) with Kj, ~ —0.4, and then boosted to the lab frame according 
to the top quark boost ((3 t ) distribution. The /3 t distribution was extracted from the 
Monte Carlo event generator directly, but can be fitted using a Gaussian distribution, 
yielding identical results (see appendix |E.2| ). Similarly, the u variable distribution 



was calculated first by distributing via a Monte Carlo technique the energy and z- 
momentum of the lepton in the top centre-of-mass frame (Ei and pf respectively), 



7 This mass/coupling combination is currently being marginally excluded by results presented by 
the LHC experiments. The aim of the present study is to assess the viability of the reconstruction 
variables. 
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Figure 6: The xtop (left panel) and u (right panel) variables for left- or right-handed Z' 
bosons decaying to ui or ut, at a 14 TeV LHC, obtained from parton- level events compared 
to a semi-analytic prediction as described in the main text. 



according to the full matrix element (see appendix |D|) and then boosting to the lab 
frame using the top quark boost distribution. The variable u was then calculated by 
taking the ratio: 

E ' + M (32 ) 
E e + &pf + E b + f3 t p z b ' v ' ' 

We ignore the mixing of helicities due to finite masses since the top quark is produced 
in association with a light quark (for which mz> S> m u ). s 

3.1.2 Simulation and reconstruction 

We will assume that the new resonances have been discovered, and that their mass 
has been measured to a satisfactory accuracy (say, O (a few %)). We will also assume 
that the spin of the new resonance has been determined by measuring the angular 
distributions of the jets originating from the u and u partons in the uu decay mode. 
It is necessary, however, to outline the details of the method for reconstructing events 
of a particular topology, for which we can form the variables we have been examining 
thus far at parton level. 

We focus on LHC proton-proton collisions at 14 TeV, in which a Z' is exchanged, 
producing a u (or w)-quark and an anti-top (or top), with a subsequent semi-leptonic 
(restricted to e or /i) decay of the top. The topology is shown in Fig. [7|. The ut and 
ut decay modes account for slightly less than ~ 50% of the total decay widths, if 
only one helicity (left- or right-handed) is present: Vz>,m=i.5 leviut/ut) = 236 GeV. 



8 These effects are small in the case of the Z' model but have been calculated in the leptoquark 
case in the following section, where the reconstruction is explained in further detail. 
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Figure 7: Production of ut from the exchange of a Z 1 and subsequent decay of the top 
into a leptonic W and a 6-jet. This mode can be fully reconstructed if one applies the W 
mass shell condition, and chooses the solution which yields the 'best' top mass. 

The leading-order cross sections for the specific topology at a 14 TeV LHC are 
a(Z' — > tu/tu — > M^u + jet) = 6.3 pb, again, if only one helicity is present, and 
including the branching ratio of the top quark to electrons or muons. We checked 
that the obvious irreducible background, SM QCD production of t + q, could be 
reduced to less than ~ 1 pb, with an appropriate invariant mass cut on the visible 
decay products of the tq pair, e.g. of around 500 GeV, with a smaller effect on the 
Z' signal, ~ 60 — 70%. We expect that other, reducible, SM backgrounds should be 
easier to reject in the helicity analysis. 

If a 6-jet is tagged and a high-p T lepton is found, along with the high-p T jet 
originating from the u, the information left missing to fully reconstruct the final 
state are the three spatial momentum components of the neutrino (assuming massless 
neutrino). We can obtain the transverse components of the neutrino momentum to a 
reasonable accuracy (C( few %)) by assuming that they are equal to the components 
of the missing transverse momentum. Then, the only remaining missing information 
is the ^-component of the neutrino momentum. In a hadron collider we do not 
possess any information on the initial ^-momentum of the system. However if we 
assume that the neutrino and the lepton originated from the decay of an on-shell W 
boson, we may apply the following mass-shell condition on their four-momenta: 



where pi and p v are the lepton and neutrino four-momenta respectively, and my/ is 
the on-shell W boson mass. This approximation is good, since the width of the W 
boson is small compared to its mass (IV — 2.14 GeV versus mw — 80.40 GeV), and 
leaves us with a quadratic equation for the ^-component of the neutrino momentum. 
To pick one of the two solutions, we choose the one that also yields a top mass closest 




(3.3) 



10 



to the on-shell top mass, via: 



m 



Ia/B = (Pl+Pu,A/B+Pb) 2 



(3.4) 



where pt, is the fe-jet four-momentum and m t ,A/B is the top quark mass obtained by 
using the solutions for the neutrino momentum p u> a/b- Once the 'best' solution is 
chosen, we possess all information required concerning the event, and we can thus 
calculate all the variables we have been examining at parton level. 

The Z' model has been implemented and has been simulated using the HERWIG++ 
event generator with initial- and final-state radiation turned on, as well as hadroniza- 
tion and the underlying event. We simulated 10 fb _1 of data, a reasonable amount 
in a near- to mid-term LHC run at 14 TeV. The events were then processed through 
the Delphes detector simulation [21], where the following minimum cuts are applied 
to the reconstructed objects: 

• PT.min for jets of 20.0 GeV. 

• Pr.min for electrons and muons of 10.0 GeV. 

The default Delphes 6-quark flat-pr tagging efficiency was replaced by a more realistic 
function of jet transverse momentum, p?, which has the form 



See for example Ref. [22] for further details. The dependence on jet pseudo-rapidity 
remained flat. An additional cut requiring the total missing transverse energy to be 
greater than 20 GeV was applied. The rest represent the default settings present in 
the Delphes ATLAS detector card, where the triggering simulation has been turned 
off. In this analysis, and the rest of this paper, we use the anti-A/p clustering algorithm 
with a radius parameter R = 0.4 to construct jets. 9 

The 'best' reconstructed top mass, after detector effects, was found to possess a 
peak at the correct on-shell top mass, ~ 174 GeV, with an approximately Gaussian 
distribution of width ~ 20 GeV. To illustrate the effect of the Delphes detector 
simulation, we show a comparison of results before and after Delphes processing in 
Figs. ^| and ^| for x top and u respectively. The comparisons between the left- and right- 
handed variable distributions after detector effects are shown in Fig. [K| It is clear 
that, at least with the minimum cuts, the u variable performs well in discriminating 
between the pure left- or right-handed top quarks. The x top variable is less different 
between the two top helicities, but the difference is still statistically significant. The 

9 It is advantageous to use a smaller radius parameter for the anti-fey clustering algorithm than 
the Delphes default one of R = 0.7, since underlying event and pile- up contaminations are expected 
to be approximately proportional to R 2 . See Ref. [23] for further details on jet algorithms and the 
underlying event. 



P(p Tjj ) = 0.08 + 0.006 x p T;s x exp (-3 x lO" 5 ^) . 



(3.5) 
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Cut set 


•^top 


u 


Min. 


40.7 


145.8 


A 


36.7 


112.2 


B 


37.6 


75.8 



Table 1: The value of X 2 /Nd.o.f. between the left- and right-handed distributions in the Z' 
model with decays to tu and tu, for the three different sets of cuts. It is evident that the 
distributions are distinguishable even for the higher cuts, with the u variable distribution 
performing best in all cases. 



advantage possessed by the variable u is clear: since no explicit reconstruction is 
required, the approximations that are associated with this process do not have a 
significant effect. As can be seen in the before/after plots in Figs. |8] and [5[ the 
u variable is less sensitive to the experimental effects, which seem to 'squeeze' the 
distributions towards middle values of the energy fraction more dramatically in the 
case of x top - 

To investigate the effect of higher transverse momentum cuts on the objects used 
in calculating these variables, we constructed two further sets of plots with higher 
cuts which we call 'A' and 'B' and are, respectively, px > 30 GeV and pr > 50 GeV 



for both jets and leptons. The resulting distributions are shown in Figs. |TT] and |12| 
where the set of cuts A is shown on the left, and set B on the right. Table |l| 
shows the value of xV^d.o.f. between the left- and right-handed distributions for the 
case of minimal cuts, as well as cuts 'A' and 'B'. 10 The value of x 2 7-^d.o.f. indicates 
how distinguishable the two distributions are statistically. It is evident that even 
in the case of higher cuts, discrimination in the specific scenario between the left- 
and right-handed modes is still possible. The higher cuts would also be beneficial 
for the rejection of proton-proton pile-up, 11 which would be detrimental at high 
instantaneous luminosity. 

3.2 Z' ->■ T+T- 

Another interesting example, which we can use to examine the effectiveness of the 
energy ratio x T , is a model of a heavy Z' possessing a t + t~ decay mode. For this 



10 For the comparison of two binned data sets, x 2 can be defined as [24] 



Ri + Si 

where Ri and Si are the number of events in i th bin of the first and second data sets respectively. 

11 Pile-up is contamination originating from multiple secondary proton-proton collisions in the 
same bunch-crossing. 
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Figure 8: Shown in the fig ures is a comparison between the results obtained for the Xtop 
variable for the 1.5 TeV flavour-changing Z' model before detector simulation (but applying 
all cuts and using equivalent jet-finding) and after the Delphes simulation for the left- and 
right- handed fermions (blue and red respectively). 
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Figure 9: Shown in the figures is a comparison between the results obtained for the u 
variable for the 1.5 TeV flavour-changing Z' model before detector simulation (but applying 
all cuts and using equivalent jet-finding) and after the Delphes simulation for the left- and 
right- handed fermions (blue and red respectively). 



model, we define a Lagrangian density similar to the one given in the previous section: 

C = g T z T , R Z'»f Rlll T R + g^ L Z'^f L ^T L 

+ g u £ R Z'»u Rlll u R + g u £ L Z'»u Llll u L + h.c. , (3.6) 

We examine discrimination between two cases, Z' — > t^t^ (g^ R = 1, g& L = 0) and 
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Figure 10: The Xtop (left panel) and u (right panel) variables for left- or right-handed Z' 
bosons decaying to ui or ut, obtained from the reconstructed events for an LHC run at 14 
TeV, with 10 fb _1 . 



Z' — > t£t l (g T zi R = 0, g^ L = 1), using the x T variable denned in Eq. (|2.2| ). 

We checked that the irreducible backgrounds to the t + t~ decay mode, e.g. orig- 
inating from the Z /'-/ — > t + t~, can be rejected using an invariant mass cut on the 
r-jet pair. For example, an invariant mass cut of ~ 500 GeV offers a rejection factor 
of 0(1O~ 4 ) for these backgrounds, bringing them to cross sections 0(0.1 pb), while 
reducing the signal only by a factor of ~ 0.7 for the case of a 1.5 TeV Z', maintaining 
a cross section of a few picobarns. Possible backgrounds to the helicity determination 
may also arise from fake r-tagged QCD events, and these would need to be assessed 
by the individual experimental collaborations. 

To compute the x T variable, one needs to reconstruct the neutrino energies from 
the r decays. For this purpose, the authors of Ref. [7] have used the collinear ap- 
proximation to reconstruct the neutrino momenta. In the collinear approximation, 
the neutrinos are assumed to be collimated to the associated r-jets. This assump- 
tion is almost always good in cases of a heavy resonance decaying to r leptons. 
Once the neutrino momentum directions are determined, the neutrino energies can 
be calculated using the missing transverse momentum constraint: 

( PmisA = ( sin 0jetl COS 0jetl SU1 fl jet2 COS jet2 \ / E vl \ 
VPmiss/ \ sin #jetl Sin0 je ti SU1 je t2 SUl 0j e t2/ \E V2 ) 

where 0j et i, 0j et i are the polar and azimuthal angles respectively, related to jet i, p^L 
are the missing transverse momentum components and E ui is the energy of neutrino 
i. 
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Figure 11: The xtop variable for a left- or right-handed Z' bosons decaying to ut or ut, 
obtained from the reconstructed events for an LHC run at 14 TeV, with 10 fb _1 , with the 
set of cuts A (left) and B (right), as explained in the text. 
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Figure 12: The u variable for a left- or right-handed Z' bosons decaying to ut or ut, 
obtained from the reconstructed events for an LHC run at 14 TeV, with 10 fb -1 , with the 
set of cuts A (left) and B (right), as explained in the text. 



However, when the two jets are back-to-back, i.e. 0j et i = 0j Ct 2 + 7T, the inverse 
of the matrix in Eq. (|3.7|) becomes singular and any small mismeasurement on the 
missing transverse energy or jet momentum directions would cause a very large error 
on the reconstructed neutrino energy [25]. The back-to-back configuration is strongly 
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Figure 13: The real part of the complex solutions that arise after applying the two 
linear missing-p^- constraints and the two quadratic r mass-shell constraints. A clear peak 
structure is seen at the input mass of 1.5 TeV. 



preferable if a heavy resonance, such as the standard model Higgs boson or a Z', is 
considered. 

One can avoid the use of the collinearity assumption by instead using information 
on the t decay vertices [26]. The most useful and best- measured attribute of these 
is their impact parameter. The impact parameter b is the displacement of a decay 
vertex in a direction perpendicular to that of the visible decay momentum, in this 
case the r jet momentum p-. Then the invisible momentum p„ must lie in the 
(b, pj) plane, so we can write p u = xh + ypj. For hadronic r decays, the invisible 
momenta are carried by single neutrinos and so their four-momenta are fixed by x 
and y for each decay. In this section we focus on hadronic r decays, by including 
a lepton veto in our event selection criteria. These four quantities are subject to 
two linear missing-pr constraints and two quadratic r mass-shell constraints, giving 
four (complex) solutions for the neutrino momenta which allow us to compute the 
invariant mass of the r pair. In Fig. O, we plot the real part of the invariant mass, 
where we use the true jets and missing transverse momenta for 10 4 events. A distinct 
peak structure is seen at the input Z' mass of 1.5 TeV. 

Unfortunately, our reconstruction method is still sensitive to the momentum 
mismeasurement. In Fig. |14], we show the hadronization and detector effects on the 
Z' mass reconstruction. In the left plot, the parton-level jets are replaced with the 
detector-level jets obtained from the Delphes simulation, in the middle plot only the 
parton-level missing transverse momentum is replaced with the detector-level one 
and in the right plot, we use all the detector-level objects. In all of the plots in 
Fig. [14], the true impact parameter was used. As can be seen, the peak structure is 
completely lost if one uses the detector objects. The event selection cuts here, and in 
all the results that follow in this section, are the default Delphes cuts for the ATLAS 
detector (pT,e > 10 GeV, prjet > 20 GeV) along with the addition of the lepton veto 
(in this case rejecting all events with an identified lepton). 

However it is important to realise that, in the model we are considering, the 
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Figure 14: The distributions of the real part of the four complex solutions resulting from 
the reconstruction of the events after using detector-level jets instead of parton-level jets 
(left plot), detector- level missing transverse momentum instead of the true one (middle 
plot) and using completely detector-level objects (i.e. both for jets and MET, right plot). 
In all plots non-smeared r vertices have been used. 

Z' mass can be measured independently from the decay modes containing the light 
quarks. In this case, we can use the mass shell constraint, (p T1 + p T2 ) 2 = m 2 z ,, as 
an extra condition to correct for the detector resolution. To accomplish this we 
construct the following likelihood function: 



J tot 



^detector 



x L 



phys) 



(3. 



where we have defined 



miss v prrnss 
E <b > 



-^detector = II - =1 (PgvP^P^) X Pj 

iph y = Pm z , x e(Re[E Ul }) x e(Re[P, 2 ]) 



(3.9) 



The p Ei , p 6i , p^ s , n 



are Gaussian probability functions centred at the origin, 



PTmiss / ' 

j.obs"i 



with arguments {E 3 -E^)/E$\ (PT m iss-p? miss )/^ b miss , (<t> PTmis - ^ 

respectively. We also use the flat probability distribution, P^ , for the (4 
with a range [-0.9:0.9] so that the probability function matches the actual probability 
of mismeasurement, which is simulated by Delphes. The likelihood L p ^ ys allows us 
to correct the mismeasured observable by requiring physical conditions with some 
probability. P mz , is a Gaussian probability function with an argument of [ra TT — rriz>), 
where mz> is the true Z' mass assumed to be measured through some other decay 
mode. B(x) is 1 if x > 0, otherwise. The probability functions we use and the ones 
that appear in Delphes are shown in appendix ||. 

For each event, we generate 1000 pseudo-events, in which the observed momenta 
are slightly shifted in a random 'direction' according to the same probability function. 
We only keep the pseudo-event, z max (corresponding to L™ x ), that provides the 
maximum likelihood. We show the rr invariant mass distribution obtained from 
Vax sample in Fig. IT3. 
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Figure 15: The rr invariant mass of the maximum likelihood sample as described in the 
text. The distribution is sharply peaked at the input Z' mass of 1.5 TeV. 



In the left panel in Fig. [16], we show the relative difference between the true 
neutrino energy and the reconstructed neutrino energy by the likelihood method. 
As can be seen, the true neutrino energy is well reconstructed on an event-by-event 
basis with about 50% error. The right panel in Fig. [16] shows the reconstructed 
x T variable in this method, using 1 fb _1 of integrated luminosity at a 14 TeV LHC 
(corresponding to ~ 7000 events before cuts). The red (blue) solid histogram is 
obtained from Z' — > t^t^ [Z' — > t£t^) sample. The lepton veto in this realistic 
case was applied by requiring no leptons with px > 10 GeV in < 2.4. The dashed 
histogram is the corresponding parton-level distribution of x T . It is obvious that 
the reconstructed x T has a very similar distribution to the parton-level one, and the 
difference between the left and right-handed x T distributions is visible even after 
the effects of detector resolution. The value of xV-^d.o.f was found to be ~ 13.8, 
indicating the high difference between the left- and right-handed x T histograms. 

3.3 Third-generation leptoquark pair-production 

Methods for reconstructing third-generation scalar leptoquark states in events where 
they were pair-produced have been studied in Ref. [13]. There, different mass vari- 
ables were constructed and all the possible combinations of decay modes were studied, 
including QCD and detector effects. The main background, ti production, was shown 
in Ref. [13] to produce negligible contribution in the suggested reconstruction vari- 
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Figure 16: The relative difference between the true neutrino energy and the reconstructed 
neutrino energy using the likelihood method described in the text (left panel) and the 
resulting energy fraction, x T distributions for the purely left- and purely right-handed 
cases (right panel) for 1 fb _1 of integrated luminosity at 14 TeV. 



ables in the vicinity of the leptoquark signal. We will focus here on the leptoquark 
types that can decay to a top quark and a tau lepton. 

For a list of states that can decay to the tr modes see Tables 1 and 2 in Ref. [13]. 
Instead of focusing on a specific leptoquark type, we consider a general scalar lep- 
toquark which possesses a branching fraction of 1 to a top quark and a r lepton. 
We call the leptoquarks of this type, with electromagnetic charge ±5/3, Sn and 
Srr, where the index indicates that the leptoquark will decay either to tiji, or tRjR 
respectively. In appendix |C] we consider leptoquarks of electromagnetic charge ±1/3 
decaying to the mixed combinations i R r L or i L r R (Srl and Slr respectively). To 
obtain results for other scenarios of leptoquarks that decay into this mode, one has 
to simply rescale the results to account for the appropriate cross sections. 



3.3.1 Parton-level results 

We produce Monte Carlo distributions of the variables outlined in section |2] for 
leptoquark states to compare to the predicted distributions at parton level. We do 
not present parton-level distributions for the r, as the corresponding parton-level 
results have been already extracted from the Monte Carlo event generator itself and 
appear in Fig. |. 

The Sxx (X G {R, L}) leptoquark can decay to tr modes, described by the 
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Lagrangian terms: 

gnRt c R T R S RR + g L Li c L T L S LL + h.c. . (3.10) 

In the present study we set either g RR = 0, gn = 1 (which we call purely left-handed) 
or vice versa: g RR = 1, #£Z = (which we call purely right-handed). The result 
for the highly-boosted x top distributions obtained for purely left- and right-handed 
events in HERWIG++ is shown in Fig. 17 with the appropriate analytic prediction. 
The proton-proton centre-of-mass energy was set to 140 TeV and the leptoquark 
mass was set to 20 TeV, so that the top quarks are well within the highly-boosted 
region. 




Figure 17: The Monte Carlo results for the energy fractions Xtop = £fe/£top for t^l (left) 
and T R t R (right) in the highly-boosted case. These are compared to analytical predictions 
as described in the text. 

We also present the top quark results originating from the decay of scalar lep- 
toquarks of mass 400 GeV at the LHC with pp centre-of-mass energy of 14 TeV. 
Figures |18| and |1^ show the left- and right-handed distributions of the variables Xt op 
and u respectively, produced using the HERWIG++ Monte Carlo event generator and 
including semi-analytical predictions. These include the effects of the finite top and 
tau masses, which introduce a mixture of helicities even though the Lagrangian terms 
are purely chiral, and the effect of the variation of the top boost in the lab frame, (3f 

The semi-analytical predictions have been produced using Monte Carlo tech- 
niques as described in section |3J]: the x top distribution was produced by distributing 
cos 9b according to P(cos9b) = (1 + P t k b cos9b), and the (3 t distribution for the top 
quark boost was extracted from the Monte Carlo event generator (see appendix |E.2| ) . 
Using either a fit or the extracted distribution yields indistinguishable results. 

Since the boost of the parent leptoquark in the lab frame and the daughter top 
quark in the lab frame are correlated, we extracted the two-dimensional distribution 
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Figure 18: The results for the energy fractions x top = £b/£top for Tl^l (left) and T^tR 
(right) for a 400 GeV mass at 14 TeV pp centre-of-mass energy. These are compared to 
semi-analytical predictions as described in the text. 
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Figure 19: The Monte Carlo results for the energy fractions u = g-qfg- for t^l (left) 
and Tutu (right) for a 400 GeV mass at 14 TeV pp COM energy. These are compared to 
semi-analytical predictions as described in the text. 



P((3 p ,(3 t ) from the HERWIG++ event generator (appendix |E.2| ). In principle, the 
f3 p distribution can be calculated for any process using the hadron parton density 
functions and an assumption for the hard process, in this case scalar SU(3) C triplet 
pair-production. The variable (3 p in the case of pair-produced leptoquarks of mass 



Mlq is then related to the centre-of-mass energy Q by (3 P = J 1 — AM^q/Q 2 . The 
method described in Ref. [4] was then used to calculate the detected polarisation of 
the top for each event. The effect reduces the polarisation on average by less than 
10%. The mixing of helicities due to the mass of the top quark is low (less than 0.1%) 
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due to the fact that it is produced along with a very light fermion (the r lepton) in 
the case we are considering. 

The u distribution was produced in a similar way, using the full polarised top 
matrix element (see appendix [D]). The W decay to a lepton and a neutrino was 
set up in the W rest frame using a polar angle 9 and an azimuthal angle for the 
lepton and neutrino momenta. These were then boosted to the top frame, where 
the b quark and W boson momenta were distributed in the top frame using a single 
polar angle 6. 12 The W mass was distributed according to a Breit-Wigner, centred 
about mw The effect of the W width was found to be small. The distribution was 
then calculated by taking the ratio: 

_ E i + Agf (3U) 

E t + /3 t pf + E b + f3 t p z b ' 1 ' ' 

where E% and Ej, are the lepton and 6-quark energies in the top rest frame and fit 
is again the boost of the top, sampled from either from the fit (Eq. (E.21|) ) or the 



Monte Carlo distribution directly. The calculations of the effect of the finite top and 
r masses in the decay of a scalar and the relation of the top axis of polarisation and 
direction of motion follow those which appear in Ref. [4] and are described briefly in 
appendix |E|. 

3.3.2 Simulation and reconstruction 

A mass reconstruction strategy for the {tr){tr) decay mode is described in Ref. [13]. 
The reconstruction there focuses on the modes S(S) — > bj(j)jiUi, S(S) — >■ b£u 3 j 2 u 2 . 
We call this the hadronic/semi-leptonic mode, as opposed to the fully hadronic mode 
which we will examine below. An important assumption, that we have already 
discussed in section |37|, which allowed for the full reconstruction of this decay mode, 
is the collinearity of the decay products of the tau leptons, owing to the fact that they 
are highly-boosted in the lab frame. This has been tested for different leptoquark 
masses in Ref. [13]. The assumption can be applied by the relation p n = ZiP^, where 
z = l,2 and the energy ratios imply that Zi > l. 13 In Ref. [13] a quartic equation was 
obtained for the energy ratio z 2 and each solution is an unique reconstruction of the 
whole event. This method provides a clean way to discriminate the leptoquark signal 
from the background. However, our preliminary attempts to reconstruct the helicity 
of this mode have indicated that both the number of events, and the quality of the 
individual four-momenta reconstruction are insufficient for detailed determination of 
the top or r helicities. Hence, we focus instead on the topology that contains two 



fully hadronic tops, shown in Fig. 20 



12 Initially there arc 9 degrees of freedom coming from the momenta of the b, i and v. Four- 
momentum conservation offers four constraints and the mass-shell conditions for the top and W 
offer a further two. This leaves us with the three degrees of freedom: 8, <fi and 8. 

13 Since the events we are considering are not a rr resonance, there is no issue with back-to-back 



r leptons as the one which previously appeared in section 3.2 
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The signal was generated using the HERWIG++ event generator, including initial- 
and final-state radiation (ISR and FSR), hadronization effects and the underlying 
event (multiple parton interactions) and the detector response was simulated using 
the Delphes package with the default ATLAS settings, modified by the 6-tagging 
function of Eq. ( |3.5| ), without the trigger simulation. Here, we also used the anti-/cr 
algorithm with radius parameter R = 0.4. The following cuts were applied on data 
corresponding to an integrated luminosity of 100 fb" 1 at 14 TeV: 

• A minimum of 6 jets (since the jets originating from the W could be identified 
as one jet in the cases where the W is highly boosted). 

• The missing transverse momentum in the event, $ T > 20 GeV. 

• Two r-tagged jets and two 6-tagged jets, all with the extra requirement that 
they have px > 20 GeV. 

Note that these are the minimal cuts that one could impose in principle experimen- 
tally, and are lower than those imposed in the reconstruction techniques proposed in 
Ref. [13]. 



Full reconstruction of the decay topology shown in Fig. |20] is possible, since we 
would only be missing the neutrino momenta originating from the decays of the 
highly-boosted r leptons, once the hadronic tops have been reconstructed. Using the 
approximation for the tau collinearity, one is left with only two unknowns, the z\ 
and z 2 energy fractions. By assuming that the x and y components of the missing 
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momentum are equal to those of the sum of the two neutrino components, we obtain 
two linear equations: 



PL SS = J#i(*i -1) + p y j2 (z 2 -l), (3.12) 
which can be solved to give: 

V T T V 

Tl 'Y) Y) TY 

^'2 "miss "j2l J miss 

Z\ = i -\ - y y - , 

PjiPj 2 ~ I'al'yi 

( YY^ f YY^ TY^ TY^ 

-, y jit miss fjlfmiss / n \ 

Z 2 = 1 - - y y ~ • 3.13 

PjlPjZ - Pj 2 Pj 2 

The invariant mass of each of the two lepto quarks may be written as m 2 s = {pt+Pr) 2 , 
resulting in the following expression: 

m 2 s = 2zip ti ■ pji + m 2 op , (3.14) 

where we have neglected the r mass term. Using Eqs. ( |3.13| ), we obtain two values of 
ms per event. Since this analysis would be performed after potential discovery, we 
would already have a measurement of the mass of the leptoquark. This would allow 
for elimination of backgrounds that may contribute and alter the energy fraction 
distributions. 

To assess the possibility of measuring the helicity of the top quarks and tau 
leptons, we generated 100 fb _1 of a fully hadronic sample for purely left-handed or 
right-handed couplings and passed them through the Delphes simulation. We then 
analysed events which contained 2 r-tagged jets and 2 6-tagged jets. We looked for 
1 or 2 jets which reconstructed the top mass in conjunction with the tops, within an 
80 GeV window. 14 For completeness, we show in Fig. I2TI the resulting reconstructed 



masses using the above method, without any attempt to optimise the resulting values. 
In the figure, we plot all the 4 entries per event, which are a consequence of the two 
leptoquark masses (ms, m§) and the two-fold combinatoric ambiguity, arising from 
the possible pairings of the top quarks and the r leptons. To obtain the best value of 
the variables x top and x T , we choose the combination which yields the best leptoquark 
mass, that is, the one closest to the true mass. We then obtain two values of the x T 
variable: x T< i/ 2 = 1/ Z\/ 2 and two values of the x top variable by using the energies of 
the two 6-jets and the reconstructed top energy. A comparison between the results 
obtained before detector simulation, but applying all cuts and using equivalent jet- 



finding, and after the Delphes simulation, is shown in Figs. ^2] and [23| on the left- and 



14 To further improve the 'top-tagging' capabilities of the analysis, one can employ a more ad- 
vanced tagging algorithm such as the one presented in Ref. [27]. For our purposes, the simpler 
reconstruction method of requiring combinations of jets to satisfy the top mass is sufficient to 
provide good results. 
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Figure 21: Mass reconstruction in the fully hadronic mode for the purely left- and right- 
handed cases for an integrated luminosity of 100 fb . There are 4 entries per event 
corresponding. 



right- handed fermions, for the variables x top and x T . The results without detector 
simulation have been normalised to the number of events resulting after Delphes 
simulation. The differences that arise at low and high energy fractions and can 
be attributed primarily to the efficiency of the tagging algorithms and the overall 
differences to the smearing of the four-momenta due to the simulation of the response 
of the detector. Comparisons of the Delphes results are shown in Fig. 23 for the purely 
left- and right-handed cases. 

To assess the discrimination capabilities of the distributions, we calculated the 
X 2 /Nd.o.f. between the left- and right-handed distributions. To investigate the effect 
of higher transverse momentum cuts, we re-ran the analysis with two higher p? 
cuts on the jets and the missing transverse momentum: set 'A-prime' with px and 
$ T > 25 GeV and set 'A' with px and $ T > 30 GeV. The results are shown in 
Table 0. 



4. Conclusions 

We have investigated variables that have been defined for the purpose of determin- 
ing the helicity of top quarks and tau leptons, in a more realistic setting than what 
has been done thus far in the literature. We first examined these analytically in 
the highly-boosted case, where the calculation is simplified, and no explicit event 
reconstruction is required for the u variable. Subsequently we focused on two spe- 
cific models: one containing a new heavy vector boson, Z', with decays to either a 
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Figure 22: Shown in the figures is a comparison between the results obtained for the 
xtop variable for 400 GeV leptoquarks before detector simulation (but applying all cuts 
and using equivalent jet-finding) and after the Delphes simulation for the left- and right- 
handed fermions (blue and red respectively). 
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Figure 23: Shown in the figures is a comparison between the results obtained for the 
x T variable for 400 GeV leptoquarks before detector simulation (but applying all cuts and 
using equivalent jet-finding) and after the Delphes simulation for the left- and right- handed 
fermions (blue and red respectively). 



light jet and a top quark or two taus, and a specific scenario in scalar leptoquark 
pair-production, in which the decay of both leptoquarks is into a top quark and a r 
lepton. We examined the flavour-changing Z' model at parton level, producing the 
relevant distributions semi-analytically and comparing these directly to the Monte 
Carlo-generated distributions. We considered experimental and reconstruction ef- 
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Figure 24: The Xtop (left panel) and the x T (right panel) variables for left- or right-handed 
tr for an integrated luminosity of 100 fb _1 modes originating from the reconstructed events 
in the fully hadronie case, using the method described in the text. 



Cut set 


•^top 




Min. 


4.0 


4.5 


A-prime 


2.9 


3.8 


A 


1.9 


2.8 



Table 2: The value of xV-^d.o.f. between the left- and right-handed distributions in the 
leptoquark tr decay scenario, for the three different sets of cuts. It is evident that the 
distributions become more difficult to distinguish for the higher cuts at the given integrated 
luminosity of 100 fb _1 . 



fects for the case of a LHC at proton-proton centre-of-mass energy of 14 TeV and 
10 fb _1 of integrated luminosity, and we investigated the applicability of the helicity 
discrimination variables. We found that in the case of a 1.5 TeV Z' that can decay 
into an up quark and a top quark, the top quark helicity can be determined, even for 
higher momentum cuts. For the Z' model with decays to t + t~ we used the r decay 
vertex information along with a likelihood method to correct for detector resolution 
effects, resulting in good discrimination between the left- and right-handed modes 
for 1 fb" 1 of integrated luminosity. For the leptoquark pair-production model, for 
400 GeV leptoquarks decaying into a top quark and a r lepton each, we examined 
reconstruction of the case when both the decaying fully hadronicaliy. Discrimination 
in this scenario is more challenging, but values of \ 2 7-^d.o.f. ~ 2 — 3 can be obtained 
even with higher than minimal cuts for an integrated luminosity of 100 fb _1 . 

To summarise, we have assessed the magnitude of the effects of QCD, cuts on 
the transverse momentum, detector effects and finally the reconstruction issues that 
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arise in this phenomenological study of helicity variables. This work is indicative of 
the difficulties that arise in 'measuring' the helicities to determine the form of the 
interactions of new particles to quarks and leptons of the third generation. To fully 
determine the potential performance of these variables, the next step would be for 
the experimental collaborations to perform similar analyses, with realistic detector 
response and particle identification techniques. 
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A. Angular variables 
A.l Definitions 

We define here two 'angular' variables, the first of which is the angle 9^, defined 
between the 6-quark and the lepton, £, in semi-leptonic top decays as shown in Fig. E5[ 
The angle is defined in the lab frame, but is shown in the figure in the centre-of-mass 
frame of the decaying top quark for illustration purposes. We consider a function of 
this variable defined by: 15 

/(cos 06,*) = 0.25 x (1 + cosfl M ) 2 . (A.l) 

We will also be considering the distance between the lepton and neutrino for semi- 
leptonic top decays, given by AR(£, v) = J~&rf[ v + 54>1 u i where 5rj£ jV and 5<j>i )V are 
the distances in the pseudo-rapidity, \t]£ — t] u \, and transverse plane angle, \4>£ — (f> u \. 

A. 2 Angular variables in the Z' flavour-changing model 

Figure |26| shows the angular variables at parton-level for the Z' model described by 
the Lagrangian density of Eq. ( |3.1|) , for a mass of 1.5 TeV at a 14 TeV LHC. In Fig. |27] 
we show the reconstructed distributions after Delphes simulation for the minimal set 



of cuts, and Figs. pj| and 29 show the corresponding reconstructed distributions for 



15 This form of the function is chosen to resemble the factor 0.25 x (1 + cos#j^) 2 , where 6£ e is 
defined in the W rest frame, which appears in the top quark decay differential cross section. 
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Figure 25: The angle Of,^, between the b quark and lepton in top decay, shown in the 
centre-of-mass frame of the top for illustration purposes. The angle is calculated in the lab 
frame throughout this paper. 



the set of cuts A and B, defined in section |3.1.2| . Table |] shows the corresponding 



X 2 /Nd.o.f. between the left- and right-handed distributions corresponding to the fig- 
ures. It is evident that these angular variables can provide equivalent magnitudes 
of discrimination between left- and right-handed top quarks as the energy function 
variables that have been used throughout the main part of the paper. 
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Figure 26: The /(cos #6^) (left panel) and A.R(£,v) (right panel) variables for left- or 
right-handed Z' bosons decaying to ui or ut, obtained from parton-level events 



B. Detector- level smearing 

The probability density functions used for the likelihood method that corrects for 
likelihood methods are shown in Fig. and Fig. |31| with dashed histograms. The 
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Figure 27: The /(cos#b^) (left panel) and AR(£,u) (right panel) variables for a left- or 
right-handed Z' bosons decaying to ui or ut, obtained from the reconstructed events for 
an LHC run at 14 TeV, with 10 ftr 1 . 
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Figure 28: The AR(£, v) variable for a left- or right-handed Z' bosons decaying to ut or 
ut, obtained from the reconstructed events for an LHC run at 14 TeV, with 10 ftr 1 , with 
the set of cuts A (left) and B (right), as explained in the text. 



solid histograms show the actual discrepancy between parton-level and detector- 
level objects as obtained by Delphes. In practice these distributions would be known 
features of the detector performance. 
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Figure 29: The /(cos 6^) variable for a left- or right-handed Z' bosons decaying to ut or 
ut, obtained from the reconstructed events for an LHC run at 14 TeV, with 10 fbr 1 , with 
the set of cuts A (left) and B (right), as explained in the text. 



Variable 


AR{£,v) 


/(cos0(M)) 


Min. 


45.7 


29.1 


A 


35.4 


19.1 


B 


34.0 


12.6 



Table 3: The value of xV^d.o.f. between the left- and right-handed distributions in the 
Z' model with decays to ut or ut, for the three different sets of cuts. It is evident that the 
distributions are distinguishable even for the higher cuts. 



C. Alternative leptoquark decay: Sxy ^bTl^lJr 

In section [0]we considered the decay of a leptoquark which we called Sxx, where 
X G {L,R}, with electromagnetic charge ±5/3 to either £.rTr or tjjri,. Here we 
consider the alternative combination of helicities, corresponding to a leptoquark with 
charge ±1/3 which we call Sxy, where 1^7 and X,Y G {L,R}. 16 The results 
for x top and x T , corresponding to the fully-hadronic leptoquark analysis constructed 
in section p\3| are shown in Fig. ^ for the minimal set of cuts. Note that in this 
case one should understand that the left-handed r should be paired up with the 
right-handed top and vice versa (i.e. red in one plot with blue in the other). 

The values of xV^d.o.f. were found to differ compared to those of the Sxx lep- 
toquark, corresponding to ~ 1.6 for the x T variable and ~ 8.2 for the x top variable. 

16 In the notation of Ref. [13], we are actually considering the decays of the upper component of 
the 51/2 doublet, s[j^. 
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Figure 30: The probability density functions for smearing of jet energy, the angle on the 
transverse plane (ft of jets, and the azimuthal angle 9 of jets. The solid histograms show 
the actual Delphes simulation results whereas the dashed one show the functions that were 



actually used for the likelihood method of section 3.2. 
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Figure 31: The probability density functions for smearing of missing transverse energy, 
the angle on the transverse plane, (ft, of the missing transverse momentum vector, and 
the actual Z' mass distribution in the model described by the Lagrangian of Eq. ( |3.6| ) 
(measured through some other decay mode than the rr mode). The solid histograms show 
the actual Delphes simulation results whereas the dashed one show the functions that were 



actually used for the likelihood method of section 3.2. 



The differences between the Sxy and Sxx leptoquarks arise due to the difference in 
acceptances of the left- and right-handed particles. These are a consequence the in- 
terplay of the efficiencies associated with the 6-tagging, r-tagging and reconstruction 
of the top quarks. 



D. Matrix element for polarised top decay 

It is useful to outline here the derivation of the polarised top decay differential width. 
The coupling of the W boson to the fermions is given by: 

9 



C 



2^2 



V//</7„(l-75)/ , W"' + h.c. 



(D.l) 
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Figure 32: The x T (left panel) and x top (right panel) variables for left- or right-handed 
tops from the Sxy leptoquark, for an integrated luminosity of 100 fb" 1 for a 14 TeV LHC. 



where g is the weak charge and Vff is the CKM matrix element corresponding to 
the fermions / and /'. The matrix element, corresponding to the diagram shown in 
Fig. |33] is given by: 

M = 9 jV tb V ff u r b (p 2 ) lfl (l - lh )W»u s tPl uJ,{q x ) lv {l - l 5 )v r f(q 2 )W»G(q 2 ) , (D.2) 

where u and v are the positive and negative frequency spinors, is the W polari- 
sation vector and G(q 2 ) is the W propagator, 



G(q 2 



(q 2 - mly) + iT w m w 



(D.3) 



Squaring the matrix element and summing over the b quark and fermion spins, we 
obtain: 



\M\=Q Tr 



-(1 + 2s^$)(f l + m top )7 fc (l - 7 5 )#j7/*(l - 7s) 



x Tr 



t x lv{l ~ 75)^ 2 7a(1 - 7s) 



x W' x W v W* x W* k 



(D.4) 



where S = l/m top (\pi\,Eipi/\pi\) is the spin 4-vector for the top-quark and Q is 
defined as: 



Q = ^\V tb \ 2 \V ff/ \ 2 x 



1 



64 



(q 2 - m 2 w ) 2 + Y 2 w m 



(D.5) 



w 



The traces can be calculated using the FORM package [28]. The first trace in Eq. ( p.4j ), 
corresponding to a top quark with spin s = ±1/2 and the bottom quarks, is given 
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\ f'(qi,n) 

Figure 33: The Feynman diagram for polarised top decay. The parentheses indicate the 
4-momentum and spin labels respectively, for each particle. 

by: 

4 {f lV \ + V \f 2 - g» k Pl ■ p 2 - e^ k M) 
+ 8 m top s{-p%S k - S^p k + g" k p 2 ■ S - e i3 » k p\S 3 ) . (D.6) 

It is obvious that the second term in the above result vanishes if we sum over s or 
set m top — » 0. The second trace in Eq. QD.4j) , corresponding to the fermions / and 
/' is given by: 

8(^2 + Q1Q2 - 9 V \i ■ Q2 + Se^qlq^ . (D.7) 

Summing over the W polarisations introduces g kX and g 11 ", and gives, for the polarised 
top matrix element squared, 

\M\ 2 (s) = 128n(p 2 ■ q x ) [( Pl - m top (2s)S) ■ q 2 ) . (D.8) 



E. Finite mass effects on top polarisation 
E.l Production polarisation 

We adapt the description given in Ref. [4] for the case of stop decay to a top and a 
neutralino, to the decay of a scalar leptoquark to a top and another fermion. The 
corresponding Lagrangian for the decay of scalar leptoquark S to tL,R and a fermion 
/ can be written as: 

C = g L hSf + g R i R Sf + h.c. , (E.l) 

where gi and g R are the left- and right-handed couplings respectively. The axis 
of spin quantisation is taken to lie along the top direction of motion in the parent 
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leptoquark rest frame. The production amplitudes for positive and negative helicity 
top quarks depend on two functions F±: 

p _ ( E t + ^top ± \pt\)(E f + m f ± \p f \) 
± ^A(E t + m top )(E f + m f ) 

where all quantities are given in the leptoquark rest frame. The functions result from 
explicit evaluation of the matrix element for the leptoquark decay using the spinor 
wavef unctions. These are given by: 

«(* s) = , v(p, s) = ( 2 ;^Z?- ) , (E.3) 

where the spinors Xs are eigenstates of the operator 1/2(7 • s with eigenvalue s, where 
s = ±1/2 corresponding to spin-up and spin-down respectively. The spinors are 
given, for s along the z-direction, by: 

X1/2& = Q , X-1/2& = (J) • (E.4) 

For the decay of a scalar to left-handed tops, for example, the matrix element is 
proportional to v(p,s)Pi,u(k,r), where p and s are the top momentum and spin 

— # 

respectively, and k and r are the fermion momentum and spin respectively. The Pl 
operator projects the upper component of u(k,r) and hence the matrix element for 
the left-handed coupling is given by: 

M L sr oc {-2s^ x l s 2s^ X l) (f%? Xr ) ■ (E.5) 
Using the relations [29]: 

(E p + m)I - a ■ p 



Vp 7 ^ 



m) 



^(E p + m)I + a.p 

^2{E p + m) 1 ' 



we obtain the following: 



V4(B< + m«.p)(B/ + ">/) 



or, writing out the matrices explicitly: 



L 2s t (E t + m top + p z Px-Wy \ (Ef + rrif -k z -k x + ik y 



D ~ s \ p x + ipy E t + m top -p z J\ k x + ik y E f + m f + kz) 

(E.8) 
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where we have defined the denominator as D = ^A(E t + m tov )(Ef + rrif). In the rest 
frame of the scalar leptoquark, we can take the z-axis to lie along the top direction 
of motion, and hence we have p x — p y — k z — k y — and p z = —k z and the matrix 
element becomes: 



M L sr oc ~^X-s ( 



E t + m top + \p z \ \ fE f + m f + \k z \ 

E t + m t0 p - \Pz\J \ E f + m f -\k z \ 

(E.9) 

Performing the matrix multiplication, one immediately notices that the functions F± 
appear as elements of the resulting matrix: 

Mi oc -2s X l s ( F + ^ Xr , (E.10) 

If we choose s = +1/2 then: 

M^cx-(0,F.)xr, (E.ll) 

or, if we choose s = —1/2 then: 

M L _ r <x-(F + ,0) Xr . (E.12) 
Upon squaring and summing over the fermion spin, r, we obtain: 

^|^ r | 2 =l-^+-! 2 «l^ + | 2 , (E.13) 

r 

and 

-£\M L _ r \ 2 = \M L _ + \ 2 <x \g L F-\ 2 . (E.14) 

r 

Since only one component of |A^g r l contributes, we may also write: 

M L ± =M L ±T , (E.15) 

and hence: 

M L ± oc g L F T . (E.16) 

The extension to M± is trivial, for which the resulting matrix elements are thus 
given by: 

Ml^g R F ± . (E.17) 

For finite m top this gives a non-vanishing amplitude for top quarks of both helicities 
even in the limit of a purely chiral vertex. It can be shown that the polarisation 
along the production axis at a parent rest frame is given by: 

\M L + + Ml\ 2 - \M L +M R \ 2 



(Pp) 



\ M L +M R\2 + \ M L +M R\2 



{\9r [ 


2 _ 


g L \ 2 )M LQ \p t 




(\9r\ 


* + \g L \ 2 )(M hQ E t - 


™>t op ) + ^9R9Lm f m top 



Xr , 



(E.18) 
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Figure 34: The effect of the sister fermion mass mf on the purely chiral top quark 
polarisation, (-Pp) c hi ra i> 111 the decay of a scalar leptoquark of mass Mlq = 400 GeV (left 
panel), or M LQ = 1000 GeV (right panel), to tf. 



where 



\Pt 



(Ml Q + ml p - to 2 ) 2 - 4M 2 Q TO t 2 op 



2Mi 



LQ 



M 2 - to 2 + 



^= Q ' ~, (E.19) 

and Mlq is the parent leptoquark mass. It is interesting to consider the effect of the 
mass of the sister fermion produced in the decay in the purely chiral case. Setting 
either g R = or g L = defines: 

(^>ch,al = (PP) {9R or aL =0) = ± M ^Et-ml P ' (E ' 20) 

If we had rrif = m top = 0, then we would get (-Pp) chiral = ±1- For a non-zero top 
quark mass, it turns out that for a wide range of values of to/, the deviation from 



2 



the value of (-Pp) chiral = ±1 is small. Figure [34] shows the effect of the fermion 
mass for the case Mlq = 400 GeV and Mlq = 1000 GeV. It is clear that when 
to j TO top , then the effect is small, and (-Pp) chiral is close to the rrif = value. 
The effect is £>(few %) even for m f ~ 100 GeV in the M LQ = 400 GeV case. For 
Mlq = 1000 GeV, the effect of the sister particle mass is of 0(few %) even for masses 
as high as ~ 600 GeV. 17 

One may also consider effects of relativistic rotation via the Wigner angle between 
the production and detection axes of the top quark [4]. However, if the velocity of 
parent particle is small or the top is highly boosted in the parent rest frame, the 
effect of the Wigner rotation is negligible. Therefore we can ignore it for our Z' 



17 Note that it may be interesting to investigate whether the helicity distributions of the top can 
provide independent information on the mass of the accompanying sister particle, especially if this 
is weakly- interacting (and hence invisible). 
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Figure 35: The fit to the top quark ft distribution for 400 GeV leptoquarks is shown 
in solid black dashes. The HERWIG++ histogram extracted from parton-level Monte Carlo 
events is shown in solid black. 




examples, since we have assumed Z' being much heavier than the daughter particles. 
Even for the leptoquark example, the effect has been checked and found not to be 
significant, since most of the events have ft, <C 1 and ft ~ 1, as we will see in the 
next subsection. 



E.2 (3 distributions from Monte Carlo 

We also show the form of the two-dimensional ft, — fit distribution for 400 GeV 
leptoquarks, extracted from the HERWIG++ event generator in Fig. For the ft 
distribution (i.e. integrated over ft) in the decay of a 400 GeV leptoquark a fit can 
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be made, shown in Fig. [35} The fit has the form of a Gaussian: 

P(P t ) = aexp(-(p t -b) 2 /c) , (E.21) 

where the parameters a, b and c were given by the fit to be a ~ 3.48, b ~ 1.03, 
c ~ 0.13. The distribution integrates to ~ 1 in (3 t e (0, 1). 
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